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We report first-principles calculations, within the density functional theory, on the structural 
and energetic properties of 180° stripe domains in (PbTi03)n/(SrTi03)„ superlattices. For the 
explored periodicities (n=3 and 6) we find that the polydomain structures compete in energy with 
the monodomain phases. Our results suggest the progressive transition, as a function of n, from a 
strong to a weak electrostatic coupling regime between the SrTiOa and PbTiOa layers. Structurally, 
they display continuous rotation of polarization connecting 180° domains. A large offset between 
[100] atomic rows across the domain wall and huge strain gradients are observed. The domain wall 
energy is very isotropic, depending very weakly on the stripe orientation. 

PACS numbers: 77.55.Px, 77.80-e, 77.80. Dj, 77.84.Cg, 31.15.A- 



I. INTRODUCTION 

The growth of superlattices composed of thin layers of 
ABO3 perovskites with different physical properties has 
become one the the most promising paths to exploit the 
coupling between instabilities in order to engineer new 
functionalities in these heterostructures.'^S For a long 
time the focus was on the electrostatic coupling between 
the layers of the superlattice,'^ and the interplay with 
the epitaxial strain.^ More recently, after the discovery of 
the appearance of a polarization from the coupling of two 
rotational modes in ultra-short period PbTiOa/SrTiOa 
superlattices,'^ the interest has evolved to include also 
the interaction between ferroelectric (FE) and antifer- 
rodistortive (AFD) modes in perovskite related layered 
compounds .EHini 

At the same time, polarization domains in ferroelec- 
tric thin films are being subject of numerous investiga- 
tions due to the recent finding of their intrinsic func- 
tional properties. One of the most remarkable is the 
conductivity at domain walls within an otherwise insu- 
lating materiali ^^ l ^^ l The great activity in the research 
of polydomain phases in ferroelectric thin films has also 
demonstrated that the domain struct ure in these sys- 
tems, with rotation of the polarizatiorP^lSI gr the for- 
mation of flux-closure structures in thin films'^^"^'^' and 
vortex structures in nano-sized fcrroelectrics,^^ prob- 
ably differs significantly from the typically assumed pic- 
ture of alternating regions where the polarization points 
along opposite directions with sharp domain walls (DW) . 

In the particular case of PbTiOa/SrTiOs interfaces, 
pioneer works were devoted to the identification of 180° 
stripe domains in PbTiOa thin films grown on thick 
SrTiOs (001) substrates. The domain structures were 
characterized both in reciprocal space (strong satellites 
around PbTiOs Bragg peaks in synchrotron x-ray scat- 
tering measurements) ,'22H25l ^nd in real space (images 
recorded by atomic force microscopy) .'^ Only lately. 



the attention has turned to the study of the domain 
structures in PbTiOa/SrTiOa superlattices. Zubko and 
coworkerJ^ have recently focused on the dependence of 
their structural and dielectic properties as a function of 
the volume fraction of PbTiOs, the electrodes, and the 
applied electric fields. While the results of the super- 
lattices asymmetrically sandwiched between Nb-doped 
SrTiOa (bottom) and gold (top) electrodes were con- 
sistent with a monodomain configuration,^* those cor- 
responding to the use of symmetrically coated SrRuOa 
electrodes (both top and bottom) suggested a polydo- 
main phase with DW motion, which dynamics might be 
quite different than the conventional one Furthermore, 
through a combination of x-ray diffraction, transmission 
electron microscopy, and ultra-high resolutio n el ectron 
energy loss spectroscopy (EELS), Zubko et al.^'^^^ have 
also recently showed the progressive transition between 
two different regimes, controlled by the thickness of the 
individual SrTiOa and PbTiOa layers. In the first regime, 
present for sufficiently thin paraelectric layers, SrTiOa 
and PbTiOs are strongly electrostatically coupled: a uni- 
form monodomain polarization is adopted throughout 
the thickness of the superlattice to minimize the depo- 
larizing field. In the second regime, when the paraelec- 
tric layer thickness is increased, SrTiOa and PbTiOs are 
decoupled: the polarization is confined within the FE 
PbTiOa layers forming domains. EELS measurements 
revealed the presence of broad interfacial layers with re- 
duced tetragonality and polarization extending over 5-6 
unit cells (u.c.) into the PbTiOa layers. Strikingly, in the 
electrostatic decoupled regime, the domain structure was 
found to be coherent over several tens of superlattice rep- 
etitions. These works pointed out the suitability of this 
system for the study of domains in ultra-thin ferroelec- 
tric films, given the behavior of the ferroelectric layers as 
quasi-independent films, while the thickness of the whole 
system prevents the charge leakage when electric fields 
are applied.*^ They also constitute a good example of 
how to tune the functional properties with respect to dif- 
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ferent factors, such as the electrical boundary conditions 
or the periodicity of the superlattice. 

On the other hand, the theoretical studies of 
PbTiOs/SrTiOa superlattices have focused, so far, on 
three different aspects of the monodomain configurations: 
(i) the analysis of the structural, elec tronic a nd ferro- 
electric prop erties of both the pristin^SHMIM] g^j^^j (jjg_ 
orderecP^I^Sl (including the presence of cation intermix- 
ing or defects) interfaces, (ii) the study of the coupling 
between the AFD instabilities and the FE polarization, 
compatible with an improper ferroelectric behaviour^ or 
(iii) the influence of strain on the previous couplingP 

In this work we perform flrst-principles calculations on 
polydomain phases of (PbTi03)„/(SrTi03)n [(n | n)] su- 
perlattices, with periodicities of n =3 and 6. Our goal 
is to gain further insight on the polarization and oxy- 
gen octahedra rotation profiles. Besides, we compare the 
differences in energies between relevant phases. The in- 
fluence of the periodicity, orientation, energy of the DWs, 
and the mixed FE-AFD-strain coupling present in these 
superlatticed^ are carefully considered. We also analyze 
the strain flelds induced by the domain structure and 
their role in the inter-layer coupling. 

The rest of the paper is organized as follows: the 
method on which the simulations are based is decribed in 
Sec. [nj In Sec, III A| we compare the energy of the differ- 
ent competing phases (monodomain versus polydomain) 
to ascertain their relative stability. The atomic struc- 
ture of the domains is analyzed in Sec. |III B[ Finally, 
the polarization profiles and strain fields are discussed in 
Sec lmCl 



II. METHODS 

The simulations have been carried out within the local 
density approximation (LDA) to the density functional 
theory (DFT) using the siesta code.'^ The rest of the 
technical parameters remain the same as in Ref. [S] 

In this work we have performed simulations of {n \ n) 
superlattices, by means of a supercell approach. Two 
values of n have been considered, n — 3 and 6, aiming 
to sample superlattices within the two distinct regimes 
experimentally observed: strong (for n ^ 4) and weak 
(for n ^ 4) electrostatic coupling between the SrTiOa 
and PbTiOs layers. 

As the starting point, an ideal structure was defined 
stacking along the [001] direction n unit cells of SrTiOa 
and n unit cells of PbTiOa. The in-plane lattice con- 
stant was fixed to the theoretical LDA value of SrTiOa 
(3.874 A). First, mirror symmetry planes were imposed 
at the central atomic layers of PbTiOa and SrTiOa, and 
an initial atomic relaxation was performed in order to 
find a reference paraelectric ground-state. Then, the ref- 
erence structure was replicated times along the [100] 
direction and Ny times along the [010] direction. Due to 
the periodic boundary conditions used in the simulations, 
Nx determines the domain periodicity, while Ny allows 
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FIG. 1; Differences in energies between polydo- 

main, monodomain and non-polar configurations in (3j3) 
PbTiOa/SrTiOa superlattices, as a function of the domain pe- 
riod Nx- Total energies of supercells are given per 5-atom per- 
ovskite unit cell. Circles represent the configurations where 
the AFD modes are not allowed {Ny = 1), while squares rep- 
resent configuration with condensed AFD modes {Ny = 2). 
Diamond indicates a configuration where the DW lies along 
the < 110 > direction, also allowing for the condensation of 
AFD modes. The monodomain phases have been labeled as 
in Ref. |H1 In the non-polar configuration, the AFD distortions 
have been considered. All energies are given with respect to 
the most stable monodomain configuration. 



to switch on {Ny = 2) and off {Ny = 1) the AFD instabil- 
ities. Following the recipe given in Ref. [13 a percentage 
of the bulk soft mode distortion was superimposed on the 
PbTiOs layers, so the polarization points upwards in half 
of the superlattice and downwards in the other half. For 
Ny — 2 superlattices, small rotations were induced by 
hand following a a'^a^c" pattern in Glazer notation. Fi- 
nally, an extra atomic relaxation of the full heterostruc- 
ture was carried out, until the maximum value of the 
Hellman-Feynman forces and the zz stress tensor com- 
ponent fefi below 0.01 eV/A and 0.0001 eV/A^ respec- 
tively [except for the (3 | 3) superlattice with N^ — 16 
and Ny — 2 (960 atoms in the simulation box) and for 
the (6 I 6) superlattice with N^ = 12 and Ny = 1 (720 
atoms), which were relaxed down to a maximum force of 
0.05 eV/A]. 

To stablish the notation, we will call the plane parallel 
to the interface the {x, y) plane, whereas the perpendic- 
ular direction will be referred to as the z-axis. 



III. RESULTS 

A. Energetics. 

For the (3 | 3) superlattices we have performed simu- 
lations of the different competing phases in order to de- 
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termine their relative stability. The energies of the poly- 
domain, monodomain and non-polar configurations as a 
function of the domain periodicity are shown in Fig. [T] 
For these superlattices, the balance between the electro- 
static energy (which tends to reduce the domain period) , 
and the DW energy density (which tends to increase it) 
results in an optimum periodicity of the domain struc- 
ture, A, of about 12 unit cells (46.5 A) (the energy for 
Nx equal 12 and 16 might be considered as equivalent 
within the accuracy of our simulations) . 

The most stable phase found in our simulations, how- 
ever, corresponds to a monodomain structure, with the 
polarization in the PbTiOa layer pointing close to the 
perovskite unit cell diagonal (configuration labeled as 
[111] in Ref. |5] and Fig. [T]). This result is consistent 
with the upturn in the domain periodicity observed by 
Zubko and coworkers,'^ suggesting that for n < 4 the 
superlattices enter into the strong coupling regime. Nev- 
ertheless, the energy difference between the monodomain 
and the most stable polydomain configuration is very 
small (of the order of 1.6 meV/5-atom-perovskite-unit- 
cell, well below the thermal energy at room tempera- 
ture), suggesting a close competition between them for 
small values of n. A small change on any external con- 
dition (growth temperature, or how fast the system is 
cooled, etc.) might help the system to overcome potential 
energy barriers, activate transitions between them and 
could, eventually, stabilize a metastable domain struc- 
ture. This fact might explain why both, polydomain 
and monodomain samples, have been observed experi- 
ment ally.f^S 

For a given domain periodicity, the energy is systemat- 
ically lowered if the rotation of the oxygen octahedra are 
allowed, with reductions ranging between 3.6 meV per 5- 
atom perovskite unit cell (for = 6) to 2.0 meV (for 
= 12). This highlights the importance of the FE-AFD 
coupling in these heterostructures. The coupling is also 
noticeable when the pattern of the oxygen octahedra in 
the polydomain phases are compared with those of mon- 
odomain configuration. We see in Fig. [2]that at the cen- 
ter of the domains (i. e. mid-distance between two DW, 
see column of atoms embodied by a bracket in Fig. [4]) 
where the polarization is purely out-of-plane, Pz, the ro- 
tations along an in-plane axis (tiltings) essentially vanish, 
resembling the case of the [001] monodomain phase re- 
ported in Ref. '51 Regarding the rotations around the 
z-axis, the larger Pz in the PbTiOs layer in polydomain 
phases penalizes the AFD modes and the rotation angles 
are smaller. The FE-AFD coupling also affects the mag- 
nitude of the polarization, resulting in a slight reduction 
of Pz when condensation of AFD modes is allowed (see 
Table 

We also find that the effect of the DW orientation 
is small: a change in the orientation of the DW from 
< 100 > to < 110 > does not affect significantly the en- 
ergy of the superlattice. This points to a rather isotropic 
DW structure, with the energy of the domains depending 
very weakly on the stripe orientation, in good agreement 
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FIG. 2: (Color online) (a) Squematic representation of the 
center of a domain in a (3|3) PbTiOa/SrTiOs superlattice 
(see region embodied by a bracket in Fig. |4|. Atoms are rep- 
resented by balls: Sr in yellow, Ti in green, O in blue, and Pb 
in grey. In panels (b)-(d) we represent the amplitude of the 
rotations (squares) and tiltings (diamonds) of each TiOs oc- 
tahedra: (b) at the center of a domain in the polydomain con- 
figuration with Nj; = 12, (c) in the ground state monodomain 
phase (with polarization in the PbTiOs layer pointing close 
to the perovskite unit cell diagonal, see Ref. 8), and (d) in a 
monodomain phase with polarization lying along [001]. 



TABLE I: Out-of-plane polarization, , at the center of the 
domains in PbTiOs/SrTiOs superlattices with Nx = 12 u.c. 
P^^'^ (PI"'"'^) stands for the polarization at the central per- 
ovskite unit cell within the PbTiOa (SrTiOa) layer. Values in 
parenthesis correspond to the root mean square polarization, 
averaged along the [100] direction. Units in /iC/cm^. 
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with experimental results,!^ phenomenological Landau- 
Ginzburg-Devonshirc theory,!^ and model Hamiltonian^ 
simulations. 

For the (6 ] 6) superlattices only one domain periodic- 
ity was simulated due to the scaling of the system size. 
For the same reason, in this case rotations of the oxy- 
gen octahedra were not allowed in the calculations of the 
polydomain phases. In view of the results for the (3 ] 3) 
superlattice, the presence of the octahedra rotations re- 
sults in a decrease of the energies of polydomains struc- 
tures and a small reduction of the polarization at the 
center of the domains. We chose = 12 u.c. (A = 46.5 
A), close to the experimental value of A = 55 A.'22l These 
polydomain phase were found to lie 3.8 meV/5-atom- 
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FIG. 3: (Color online) Local polarization profile of polydo- 
main structures in (PbTi03)„/(SrTi03)n superlattices with 
(a) n = 3 and (b) n = 6. The PbTiOs and SrTiOs are 
depicted as gray and white regions respectively. Red dashed 
squares in the SrTiOs layers mark the position where antivor- 
tices are formed. 



perovskite-unit-cell below the most stable monodomain 
phase without AFD distortions. 

The change in the most stable phase [from mon- 
odomain in the (3 | 3) to polydomain in the (6 | 6)] 
indicates a crossover between the weak and strong elec- 
trostatic coupling regimes described in the Introduction. 
Experimentally this transition was inferred to occur grad- 
ually, with a minimum of the domains size observed at 
n ~ 4 (Ref.EIJ. In related KTaOa/KNbOa superlattices, 
the critical periodicity for the crossover ranges between 
7 < n < 15 (experiments from Ref . . and 12 < n < 24 
(shell models simulations from Ref. [37]). 



B. Atomic structure of the domains. 

Figure [3] shows the local polarization profile of the most 
stable polydomain configurations found for the (3 | 3) 
and (6 | 6) superlattices (both with a domain periodicity 
of Nx — 12 u.c). Similar patterns are obtained for other 
domain sizes. In order to be able to make quantitative 
comparisons, from now on only simulations where oxy- 
gen octahedra were not allowed are discussed. The local 
polarization is obtained calculating the polarization of a 



unit cell centered on every cation of the system (except 
at the interfaces where no "bulk-like" unit cell can be 
chosen) and using the displacement of the atoms with 
respect to the ideal phase and the bulk Born effective 
charges of the corresponding material, cither PbTiOs or 
SrTiOa, depending on the layer the cation belongs to. 

Near the DWs the local polarization pattern clearly 
displays a continuous polarization rotation within 3 u.c. 
around the DW, connecting two 180° domains. 

As suggested by the experimental results exposed in 
Ref. 1311 the examination of Fig. [3] reveals that the ac- 
tual domain structure in this kind of systems is often an 
intermediate case between 180° domains and the closure 
domains commonly found in ferromagnets, displayin g ro- 
tation of the polarization upon approaching the VW^^^ 
(although the length scale over which the polarization ro- 
tation takes place is only a few unit cells in the case of 
ferroelectrics, in contrast with the several nanometers or 
even microns typical of ferromagnets) . However, we have 
to keep in mind that in ideal closure domains, the diver- 
gence of the polarization vanishes everywhere and, there- 
fore, the depolarizing field is perfectly screened. In our 
simulations, the polarization of the SrTiOa indicates the 
presence of a residual depolarizing field and thus, strictly 
speaking, our domains do not constitute perfect domains 
of closure. 

Our results also support the robustness of the ro- 
tation of polarization and the formation of vortices 
in ferroelectric nanostructures suggested by previous 
theoretical studies. These geometries have been pre- 
dicted to exist independently o f (i) the used method- 
ology (including first-principleSj'iSESl niodel hamiltoni- 
g^jjgjfflHa pjjg^gg f^eld models,''^ and phen omenological 
Devonshire-Ginzburg-Landau theories^SES) ^ and/o r (ii) 
the electrostatic boundary co nditio ns (with metallicP^EH 
or semic onduc ting electrodes,'^^'^ or even in free stand- 
ing slabs.l22E2l) 

It is remarkable to see that the polarization rotation in 
the PbTiOa/SrTiOa superlattices is mostly due to large 
in-plane displacement of the Pb atoms at the PbO lay- 
ers in the vicinity of the interface. This contrasts with 



the predicted domains in BaTiOa/SrRuOa capacitors 



where the in-plane polarization is due to the displace- 
ments of the Sr atoms in the first layer of the electrode. 
Here, the Pb atoms move of the order of 0.2 A, a displace- 
ment large enough to be detectable with the recently de- 
veloped atomic-resolution aberration-corrected transmis- 
sion electron microscopy. Using this technique, polariza- 
tion rotation at DWs have been experimentally observed 
in ferroelec tric th in films with thicknesses of a few tens 
of unit cells .'i^'^Sl However the high quality level achieved 
during the last years in the growth of short-period super- 
lattices, together with the large in-plane displacements 
predicted, make this kind of system particularly suited 
for the observation of the formation of vortices at do- 
main walls in ultrathin films, comparable in size to the 
simulated systems listed above. 

Interestingly, within the SrTiOa layer and close to the 
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DW we do observe the formation of antivortices; a local 
polarization pattern where two dipoles point face to face 
and two tail to tail (see red dashed squares in Fig. |3] 
These antivortices have also been recently predicted to 
form in epitaxial BiFeOs films!^ 



C. Polarization profiles and strain field 

Within a polydomain configuration, there is no need 
anymore to keep constant the normal component of the 
polarization at the interface, Pz, since the electric fields 
that arise from its discontinuity are efficiently screened 
by the presence of domains. As a consequence Pz, that in 
the monodomain configuration is continuous throughout 
the superlattice, in the polydomain case is very inhomo- 
geneous with polarization mismatches at the center of 
the domains of 34 /xC/cm^ and 54 /zC/cm^ for the (3 | 3) 
and (6 | 6) superlattice respectively (see Table |l]). 

The layer-by-layer polarization of Fig. [sja) and [sjc) 
shows that, within SrTiOs, it converges to a rather homo- 
geneous well defined value, that at the center of the do- 
main decreases from 31 fiC/cm^ (for n = 3) to 21 /xC/cm^ 
(for n = 6). On the contrary, the PbTiOs layer displays 
a smooth variation of the polarization, with a progres- 
sive reduction spanning over a length of three unit cells 
into the PbTiOs layers from the interface. The great re- 
duction of the polarization of the SrTiOs layer upon a 
increase in n, together with the out-of-plane polarization 
at the center of the PbTiOa layer rapidly approaching 
the bulk value (83 ^C/cm^), again supports the gradual 
electrostatic decoupling of the ferroelectric layer.'Sii 

In Fig. [5] we also plot the variation of the local tetrag- 
onality across the superlattice, calculated for the same 
perovskite unit cell surrounding each cation. The layer- 
by-layer tetragonality averaged along the [100] direction, 
plotted as black squares in Fig. [sj^b) and (d), displays 
a variation that is well correlated with that of the po- 
larization: an almost constant value inside the SrTiOa 
layer and a smooth increase from the interfaces towards 
the center of the PbTiOa. The polarization of the SrTiOa 
layer induces a slight tetragonality of this material. This 
reduction of symmetry with respect to the cubic unit cell 
of bulk SrTiOs is consistent with the decrease of the ^23- 
Cg splitting observed in the EELS spectra of this system 
by Zubko and coworkers.*^ 

Besides, in Fig. [sj^b) and (d) we show the layer-by- 
layer tetragonality across the superlattice at the center 
of a domain with the polarization pointing up. Strik- 
ingly, and contrary to the polarization, the variation of 
the local tetragonality as we move across the superlattice 
is very asymmetric. The tetragonality reaches its maxi- 
mum value in the PbTiOa layer at the bottom interface 
with respect to the polarization direction, and gradually 
reduces its magnitude as we move towards the top inter- 
face. In the SrTiOs layer, the strain gradient is smaller 
and, forced by the imposed periodic boundary conditions, 
it has opposite sign. 
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FIG. 4: (Color online) Schematic view of the (3 j 3) superlat- 
tice with TVj, = 12 and Ny — 1 indicating how local values of 
the magnitudes plotted in Fig. [5]and[7]are defined. Red (blue) 
lines represent local values of in-plane, a, (out-of-plane, c) lat- 
tice constants, measured from the in-plane (out-of-plane) dis- 
tance between equivalent cations of the same chemical specie 
in consecutive unit cells along the x (z) direction. Magni- 
tudes with subscript 1 (2) indicate unit cells centered on a 
[001] AO (Ti02) atomic plane. Local polarization is marked 
with arrows. Black dotted lines indicate the offset between 
[100] atomic rows to the left and right of the domain walls, 
defined as the relative vertical shift of A-cations in a given 
atomic plane. Bracket at the bottom of the up domain indi- 
cates the position of its center, where the values plotted as 
empty symbols in Fig. [5] and [7| are obtained. Finally, domain 
walls are represented by dashed lines. 



The analysis of the polarization and tetragonality pro- 
files reveals that the formation of domains in the su- 
perlattices is associated with complex distortions. The 
characteristics of the strain field in this system can be 
explained as a combination of different effects. 

On the one hand, in PbTiOa the off-center displace- 
ments of both the Pb and Ti cations contribute to the 
polarization. Therefore the Pb atoms displace along z 
in opposite direction in the up and down domains. This 
gives rise to an offset between [100] atomic rows to the 
left and right of the DW [see Fig. [6]^ a)]. A sizable offset of 
0.6 A was already predicted by Meyer and Vanderbilt in 
180° stripe domains in bulk PbTiOa.^S As in Ref. 47, for 
the PbTiOa/SrTiOa superlattices we quantify this offset 
for a given layer as the difference in the z-coordinate of 
a equivalent A-cation at the center of opposite domains 
(see Fig.|4|. The layer by layer offset, shown in Fig. [Tj^a) 
and (c), amounts up to almost 0.5 (0.45) A at the middle 
of the PbTiOa layer in the (6 | 6) [(3 | 3)] superlattice. 
Although the offset of opposite domains is partially ac- 
commodated by the interfaces - which reflects in the in- 
crease (decrease) of the tetragonality at the bottom (top) 
interface in Fig. [5jb) and (d) -, it still propagates into 
the SrTiOa, amounting a sizable ~ 0.1 A. 

On the other hand, in thin films the polarization in 
PbTiOs rotates at the DW. Indeed, as pointed out above, 
our simulations show that large in-plane displacements of 
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FIG. 5: Left panels: layer-by-layer out-of-plane polarization, 
Pz, inferred from the Born effective charges and the atomic 
displacements for (a) a (3 | 3) and (c) a (6 | 6) superlattice. 
Right panels: layer-by-layer tetragonality for (b) a (3 | 3) and 
(d) a (6 I 6) superlattice. Empty symbols represent values at 
the center of an up domain, while filled symbols correspond to 
averaged values (root mean square in the case of polarization) 
along the [100] direction. 



the Pb atoms (up to 0.2 A) take place a the interfaces. 
It is sensible to argue that this in-plane polarization is 
coupled with an in-plane strain and, as it is schemat- 
ically depicted in Fig. |6jb), it pushes the DW in the 
same direction of the polarization. This effect is rein- 
forced as consecutive DWs become closer, as it happens 
in ferroelectric thin films [Fig. [6|^b)]. As a consequence, 
the in-plane lattice constant is expected to be enlarged 
at the top interface (with respect to polarization direc- 
tion) and compressed at the bottom interface. To test 
this hypothesis we have performed a detailed analysis of 
the strain field in the system, calculating for every indi- 
vidual perovskite unit cell the local values of the in-plane 
lattice constant, a (see Fig. |4] for an explanation about 
how it is computed). The local in-plane strain, calcu- 
lated as Ell = a/ao — 1, where oq — osrTiOa — 3.874 A 
is plotted in Fig. [^b) and (d) for the (3 | 3) and (6 | 6) 
superlattices respectively. It shows a variation with re- 
spect to the position along the z-direction that can be 
clearly correlated with that of the tetragonality, shown 
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FIG. 6: Schematic representation of the distortion induced 
by the domain structure in (a) bulk PbTiOa, (b) PbTiOa thin 
films and (c) PbTiOa/SrTiOa superlattices. (a) In bulk, dis- 
placements of Pb cations cause an offset between [100] atomic 
rows across the DW. (b) In thin films, in addition to the offset 
between domains, rotation of the polarization near the inter- 
face is responsible of a non- vanishing strain gradient 
(c) In the case of the PbTiOs/SrTiOs superlattices, the off- 
set and modulation of the strain field in the PbTiOs layer (in 
grey) propagates into the SrTiOs (in white). 



in Fig. ^h) and (d) : PbTiOa unit cells close to the bot- 
tom interface (with respect to the polarization direction) 
are compressed in-plane, and as a result, they tend to 
elongate along the z axis. Conversely, at the top inter- 
face, the material is expanded in-plane and presents a 
reduced tetragonality. Since the in-plane polarization is 
confined at the interfaces, large strain gradients ^^j- can 
be anticipated. In fact, huge values are obtained from 
our simulations: up to 4 • 10^ m~^, more than seven or- 
ders of magnitude larger than those obtained in bending 
experiments on SrTiOs'*^. Similar distortions are found 
in the SrTiOa layer, although with opposite sign of the 
strain gradients. 

Finally, in the superlattice, the SrTiOs atomic layers 
closer to the interface are forced to follow the local in 
plane expansion or contraction of the PbTiOa [see Fig. 
|6][c)]. However, the presence of an inhomogeneous strain 
in the SrTiOa is energetically very costly. Therefore it 
rapidly recovers a nearly constant in-plane lattice param- 
eter [see Fig. ^d)]. 

The combined effect of the domain offset and the mod- 
ulation of the strain field can explain the large out-of- 
plane coherence length of domain structure previously 
observed in this superlattice.^! Even for relatively large 
thicknesses of SrTiOs, for which the ferroeletric layers 
can be considered as being electrostatically decoupled, 
domains with the same orientation of the polarization 
are aligned along the [001] direction. 

The large strain gradients associated with the poly- 
domain configurations might, in addition, significantly 
affect the local polarization pattern through the flexo- 
electric effect .BSHISl According to the sign of the gradients 
relative to the polarization direction the flexoelectric ef- 
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FIG. 7: Left panels: local layer-by-layer offset between [fOO] 
atomic rows to the left and right of the DW for (a) a (3 | 3) and 
(c) a (6 I 6) superlattice. Right panels: local in-plane strain 
across the center of an up domain in (b) a (3 | 3) and (d) a 
(6 I 6) superlattice. A large non-diagonal component of the 
strain gradiend, ^^|^ can be observed close to the interfaces. 



feet tends to increase the p olariz ation in the PbTiOa and 
decrease it in the SrTiOs IMlss] numerical quantifica- 
tion of the corresponding non-diagonal component of the 
flexoelectric tensor is extremely subtle, since it might be 
hidden by strain contributions via piezoelectric effects, 
and is out of the scope of this paper. However, assuming 
a flexoelectric coefficient of the order of I nC/m (typical 
of ferroelectric perovskite^^ll) ^ and with the strain gradi- 
ents extracted above, we estimate that the flexoelectric- 
induced polarization can reach values of a few /iC/cm^ 
at the interfaces or at the center of the short period su- 
perlattices. 



IV. CONCLUSIONS 

In summary, using accurate first-principles simulations 
we have studied the domain structures in short-period 
(PbTi03)„/(SrTi03)„ superlattices. The most impor- 
tant conclusions that can be drawn are: (i) the domain 
structures might compete in energy with monodomain 
configurations, (ii) The domains are rather isotropic, a 
change in the orientation of the DW from < 100 > to 
< 110 > does not affect significantly the energy of the 
superlattice. (iii) From the structural point of view, they 
display polarization-rotation, similar to the one theoret- 
ically predicted in ferroelectric nanocapacitors and re- 
cently observed in various ferroelectric ultrathin films, 
(iv) Our results suggest the progressive transition as a 
function of n from a strongly electrostatic coupled regime 
(where the ground-state is a monodomain configuration 
with a constant out-of-plane component of the polariza- 
tion preserved throughout the structure), to a weakly 
coupled regime (where the polarization is confined within 
the PbTiOa layers forming domains), (v) The evolu- 
tion of the out-of-plane layer-by-layer polarization and 
tetragonality within the SrTiOs and the PbTiOs layers 
are consistent with the t2g-eg splitting inferred from the 
unit-cell-resolution recently measured by EELS exper- 
iments.'^ (vi) Large offsets between [100] atomic rows 
across the DW and huge strain gradients (seven orders 
of magnitude larger than those obtained in bending ex- 
periments on SrTiOj^S]) are observed. The contribution 
of both of them can be responsible of the out-of-plane co- 
herence of the domain structure found experimentally.'^ 

This knowledge should complement the experimental 
studies and could encourage the design of new artificial 
structures with even more appealing functionalities. 
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